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Based on a simplest molecular orbital theory of H^, a three-parameter model potential function 
is proposed to describe ground-state diatomic systems with closed-shell and/or S-type valence-shell 
constituents over a significantly wide range of internuclear distances. More than 200 weakly and 
strongly bound diatomics have been studied, including neutral and singly-charged diatomics (e.g., 
H2, Li2, LiH, Cd2, NaJ, and RbH~), long-range bound diatomics (e.g., NaAr, CdNe, He2, CaHe, 
SrHe, and BaHe), metastable molecular dications (e.g., BeH^"*", AlH"^"'', MgJ^, and LiBa^"^), and 
molecular trications (e.g., YHe++''' and ScHe+++). 
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Modeling the interaction potential of diatomic systems 
is of fundamental importance to many issues 0, I2i S 
0,0, 0,0, 13 J including atom-atom collisions, molecular 
spectroscopy, prediction of cluster structures, molecular 
dynamics simulation, chemical reactivity, matter-wave 
interferometry, and transport properties for more com- 
plex systems. Also of great interest are the potential 
functions for long-lived metastable doubly- or multiply- 
charged ions that are relevant to high-density energy 
storage materials and to characterization and analytical 
methods for biosystcms. 

Modern spectroscopy, diffraction, and scattering tech- 
niques Q i provide a direct experimental approach to 
studies of interaction potentials of diatomic systems. In 
particular, diatomic potentials can be inferred from the 
spectroscopy data by three general approaches ^10]: (i) 
the Wentzel-Kramers-Brillouin (WKB) Rydberg-Klein- 
Rees (RKR) method, (ii) the WKB-based Dunham ap- 
proach, and (iii) semiempirical or empirical procedures. 
On the theoretical side, a diatomic potential curve may 
be predicted directly by ab initio calculations and 
quantum Monte Carlo simulations These theoret- 

ical methods can, in principle, be very accurate when 
sufficient electronic configurations are included in the cal- 
culations, but can be prohibitively expensive in weakly 
bound systems and/or many-electron systems "s]. 

Numerous attempts to analytically model diatomic 
potentials have been made jj, 0, 0, 0, E, 
lla |. The well-known potential functions include Morse, 
Born-Mayer, Hulburt-Hirschfelder, Rosen-Morse, Ryd- 
berg, Poschl- Teller, Linnett, Frost-Musulin, Varshni III, 
Lippincott, Lennard-Jones, and Maitland-Smith poten- 
tials 0, HqL as well as the celebrated Tang-Toennies 
potential m and the recently proposed Morse-based 
potentials These potentials usually aim to de- 

scribe either strongly or weakly bound, neutral or singly- 
charged diatomics and often lose their validity for ei- 
ther small or relatively large internuclear distance (de- 
noted R hereafter). Thus, recent effort has been de- 
voted to the construction of hybrid potentials, which 
use different functions for different interaction regions of 



i? 0,113 El El m and thereby need more than four po- 
tential parameters. Well-known examples of hybrid po- 
tentials include the combined Morse- van der Waals [lOl| . 
general Buckingham- type exp(n,m) E3i Cvetko E3j and 
Bellert-Breckenridge 0| potentials, as well as the most 
recently proposed Rydberg-London potential |16| . For 
metastable doubly- or multiply-charged molecules, none 
of the above-mentioned potential functions is able to de- 
scribe their ground states. To date, only few theoret- 
ical models J^J that were specifically designed for 
metastable molecular dications |^ have been proposed. 

The goal of this Letter is twofold. First, we propose a 
molecular-orbital theory based approach to obtain a very 
simple analytical potential of diatomic systems. The po- 
tential function thus obtained has significant applicabil- 
ity insofar as it can describe a wide variety of diatomic 
molecules (including both weakly and strongly bound 
systems) with good accuracy for almost the whole range 
of R but excluding the large-i? limit. Second, we show 
that this potential function can also describe metastable 
doubly-charged diatomics as well as singly- and triply- 
charged ones. Specifically, we advocate a very simple 
three-parameter ground-state potential function that is 
applicable to more than 200 diatomics with closed-shell 
and/or S-type valence-shell constituents (atoms or ions 
whose shells are closed or whose valence shells are S- 
orbital). These include neutral and si ngly -charged di- 
atomics, long-range bound diatomics metastable 
molecular dications and molecular trications [2ll |. 
The details for these systems and the associated param- 
eters of our model potential are given in Ref. . 

We require a few-parameter potential function to sat- 
isfy the following basic conditions: (i) Its asymptotic 
value Eoo for i? ^ cx) is finite, (ii) A global poten- 
tial minimum i^min at the equilibrium distance Rg is al- 
lowed, (iii) It approaches infinity as i? — > 0. (iv) One 
local potential is allowed to de- 

scribe metastable systems, (v) Both Coulomb and ex- 
change interactions can be described by using only few 
parameters. To seek such a potential function we re- 
visit the molecular-orbital theory 0| as applied to HJ, 
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the simplest single-electron diatomic system, with the 

associated Hamiltonian H — xj'^ —— — + -h (in 

atomic units), where rA and tb denote electron- nucleus 
distances. This case can be solved exactly, but here it is 
used as a reference system to understand how the sim- 
plest version of the molecular-orbital theory may be im- 
proved. To that end consider the S-type trial function 
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ofH+: ^ = ci\^i)+C2 



where 



"0/ 



^ (the Is- 



orbital of H atom) . The energy of the bonding orbital is 
then given by 



E{R) = 



JiiR)+K,{R) 
l + So{R) 
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. In the literature [3, El] , Ji and Ki are called the 
Coulomb and exchange integrals, respectively, and Sq 
is the overlap integral between the orbitals |(/)^) and 
|(/)^). Figure 1(a) shows the resultant potential curve 
T+ V -0.56483 hartree, 



of Hj . The minimum energy E 
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located at Re = 2.500 bohr. This should be compared 
with the most accurate data Emin = —0.60263 

hartree, at Re = 1.999 bohr. Clearly then, while the 
analytical potential function of derived above sat- 
isfies most of the general pair potential requirements 
set above, quantitatively it should be improved. In- 
deed, if polarization and even diffuse functions are in- 
cluded in the trial function, then the potential curve 
in the bonding region has a much better performance. 
As illustrated in Fig. 1(b), by using couple cluster 
method with single and double excitation (CCSD) '2^ 
with ST0-3G (Is orbital only), 6-31G(d,p) (including 
polarization function), and 6-311-1— |-G(3df,3pd) (includ- 
ing diffuse functions) Gaussian- type basis sets, one ob- 
tains Emin = -0.582697, -0.594490, -0.602207 hartree 
at Re = 2.004280, 1.948160, 1.999899 bohr, respectively. 

We now introduce a simple analytical potential func- 
tion to improve the above potential for Hj. That is. 



E{R,a,f3,j) = E^ 



MR,j)+K,iR,a,l3) 
l + SoiR) 



(2) 



where parameter 7 is introduced in the Coulomb inte- 
gral Ji, i.e., Ji(i?, 7) = e~^^^ (1 -I- ■^), and two parame- 
ters a and (3 are introduced in the exchange integral Ki , 
i.e., Ki{R,a,l3) = e""^ (-^ - PR). Below we briefly dis- 
cuss the meanings of the three parameters in the light 
of the polarization approximation A detailed discus- 
sion of this issue is presented in Ref.|2^. In the first- 
order polarization approximation, Eq. (2) can be rewrit- 
ten as E{R,a,f3,j) = E{Ep,eex) = Ep-[l- 5o(i?)]eex, 



where Ep — E^o + Ji{R, 7) and ee 



So(R) 



'-^[Ji{R,l)- 



Ki{R,a, /3)] are the polarization and exchange energies, 
respectively (For one-electron H^, the exchange energy 
can be interpreted as resulting from the electron hop- 
ping back and forth across the median plane between 
two protons |5l| , therefore refering to the exchange of two 
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FIG. 1: The potential energy curve of the ground state of 
HJ: (a) Equation (1) (dot-dashed line) and Equation (2) 
(solid line, a = 1.0511106, /3 = 0.917034242, 7 = 2.25 
); (b) Equation (1) (dot-dashed line), CCD/STO-3G (dot- 
ted line), CCSD/6-31G(d,p) (dashed line), and CCSD/6- 
311-f-+G(3df,3pd) (solid line). The filled dots in (a) and (b) 
are the most accurate data reported in the literature l23l . The 
inset in (a) is for the short-range region. 



protons). Clearly, parameter 7 directly adjusts Ji(-R, 7) 
and hence the polarization energy Ep. Because Sex also 
depends on Ji(i?, 7), the introduction of 7 also affects 
the dispersion (positive) part of eex- Through the term 
Ki{R, a, P), parameters a and /? are used to account for 
the i?-dependence of eex that is already affected by 7. In 
particular, the induction part (the negative term) of eex 
is adjusted only by parameter a, and parameter /3 fur- 
ther adjusts the dispersion part of eex through the neg- 
ative term of Ki{R,a, P). Certainly there are alterna- 
tive approaches for realizing these adjustments, but the 
new potential function constructed above includes both 
the Pauli repulsive term and the well-known Born- 
Mayer "exponential" form Ae~^^. This is different from 
Tang-Toennis 0, Cvetko [13, and Rydberg-London [l| 
potentials, whereas only the Born-Mayer form appears 
as their repulsion terms. It should also be stressed that 
although E{R, a, /3, 7) now has three adjusting parame- 
ters, it is still analogous to Eq. (1) in many aspects (e.g., 
satisfying all the pair potential requirements set above). 
Based on this three-parameter potential function, we find 
that the potential curve for Hj, as shown in Fig. 1(a), 
would agree very well with the most accurate data avail- 
able in the literature [23| if we choose a = 1.0511106, 
P = 0.917034242 and 7 = 2.25. This confirms that a, 
P, 7 can be properly adjusted such that contributions 
of both the polarization and exchange energies can be 
accounted for in an efficient way, thereby achieving, in 
effect, the same goal as that of using larger basis sets (po- 
larization, diffuse functions) in the trial wavefunctions. 

Certainly our real motivation is to extend this sim- 
ple and successful procedure from to other multi- 
electron diatomic systems. A number of established re- 
sults about the electronic structures of diatomic systems 
suggest that this is possible for ground-state diatomics 
with closed-shell and/or S-type constituents (see details 
in Ref. 22] ). In particular, in the zeroth-order approxi- 
mation, the outermost electrons in a multielectron sys- 
tem move in the Hartree-Fock self-consistent field or the 
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FIG. 2: The comparison between the new potential (dot- 
dashed line, a = 1.5065756, /3 = 2.48475652, 7 = 1.45), 
11-parameter-fit model potential (solid line, R.ef . l26|| ) . hy- 
brid Rydberg-London potential (dashed line, Ref. |lq|'l. Morse 
potential (dotted li ne) and the most accurate ab initio data 
(filled circles, Ref. (lUl 'l for hydrogen molecule H2. Inset in 
the rightmost figure is the enlarged part between 3.0 and 5.0 
bohr. 
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FIG. 3: The potential energy curve of the ground state of 
BeH++ (a = 0.687, /3 = 1.43632004, and 7 = 0.1185) and 
A1H++ (a = 0.585984, (3 = 0.796691521, and 7 = 0.0365} 
The filled/open circles denote the previous data from Ref. 



effective potential of all the core electrons and the pos- 
itive nucleus, and the asymptotic exchange energy of a 
multielectron system can arise primarily from the outer- 
most electrons. The exchange interactions between two 
multielectron atoms, which play a crucial role in chemi- 
cal bonding, are dominated by the exchange of a single 
pair of electrons, and the associated exchange energy is 
given by that of a single electron pair multiplied by a con- 
stant. Based on the polarization approximation |^], the 
ground-state potentials E{Ep,eex) of H2 and other mul- 
tielectron diatomic systems, when expressed in terms of 
the polarization and exchange energies, can take a sim- 
ilar form to that of HJ, despite that their origins 
of the exchange energy are totally different. Motivated 
by these known theoretical results, we have carried out 
extensive studies of more than 200 diatomic systems for 
which experimental or ab initio data are available. We 
find that, indeed, the above three-parameter potential 
can be proposed as a widely applicable potential func- 
tion for ground-state diatomics with closed-shell and/or 
S-type constituents. 

To determine the three parameters of the proposed po- 
tential function we suggest several numerical approaches 
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FIG. 4: The computed vibrational levels (filled circles) of 
■^LiH ( a = 0.8885591, l3 = 1.51479003, 7 = 0.345), H2, 
Na''°Ar (q = 0.6365, /3 = 0.0229701414,7 = 0.50), and 
ii3Cd^0Ne ( a = 0.89314, /3 = 0.33243839, 7 = 0.40). The 
solid horizontal lines denote the measured vibrational levels 
for Na*°Ar :2g3 and CdNe (suggested values from Ref. ^ 
based on the average mass of CdNe). The open circles (dashed 
lines) are RKR potential points (measured vibrational levels) 
for H2 301 and ^LiH 'j^. 



in the Appendix C of our supplementary material p3|. 
The model potential curves thus determined for more 
than 200 weakly and strongly bound diatomic systems 
agree with the available experimental or theoretical data, 
with the agreement in many cases much better than 
one could naively anticipate from a three-parameter po- 
tential (see Table 1 and Fig. A in Ref.[23). Below 
we discuss some sample results. In particular. Fig. 2 
shows that the potential curve for H2 is in good agree- 
ment with the recent 11-parameter model potential 
and the most accurate data 11], thereby giving a bet- 
ter performance than the Morse function |3| an d the 
most recent hybrid Rydberg-London potential . Even 
more significantly, our potential function is applica- 
ble to metastable S-type molecular dications |2y| (e.g., 
He++, Be++, BeH++, Mg++, MgH++, BH++, A1H++, 
LiBa"*""*", KBa"'""'", NaBa"*""*", and BaJ"''), and molecular 
trications (e.g.,YHe+++, ScHe +"'""'") as well as neu- 
tral and singly-charged diatomic systems. The poten- 
tial curves for BeH++ and A1H++ using our potential 
function are shown in Fig. 3, where the potential bar- 
riers agree well with the literature dataj^^. Figure 4 
displays the calculated rotationless vibrational levels for 
^LiH, H2, CdNe and Na^^Ar (see Ref. for results 
of isotopes), reachin g go od accuracy as compared with 
experiments 125 123. l30l|. Quite unexpectedly, even 
for very weakly long-range bound diatomics such as 
^Hcz, ^^Ca^He, ^^Sr^He, and ^^'^Ba'^Re, we are able to 
find a set of potential parameters that predict a sin- 
gle vibrational level at -0.107, -67.099, -59.875, and - 



4 



48.560 /ieV consistent with the recent hterature 

data, -0.0999, -67.161, -59.573, and -48.279 fieV, respec- 
tively dip. 

For the metastable dications He^"*", Bc^^, BcH++ 
and MgJ^, we found that they can support 5, 18, 8 
and 20 vibrational levels that again agree with previous 
studies [23, ■ Furthermore, with the new potential 
function we predict that the metastable dication A1H++ 
can support 12 vibrational levels. The estimated life- 
times for the lowest four vibrational states of BeH"'""'" are 
T = 4.9 X 10^°, 3.3 X 10^ 4.8 x lO'^, and 130 ^is, and 
those for the lowest six vibrational states of A1H++ arc 
r = 2.8 X 10^6, 1.8 X lO^^, 2.0 x 10^°, 3.3 x 10^ 8.3 x lO"* 
and 288 fis (see Table 8 in Ref.01). Note that BeH++ 
and A1H++ have been recently observed to survive flight 
times of about 4 and 7 fis, respectively [33jj, thus sup- 
porting our calculations. 

Before concluding we make one final mark. In the 
large-i? limit where the atomic electron clouds do not 
overlap considerably, the interaction energy of an atomic 
pair is given by the well-known multipolar dispersion ex- 
pansion J2n=3 C2n/i?^" Illllll. In this hmit our model 
potential approaches Eoo exponentially, a feature differ- 
ent from that suggested by the multipolar dispersion ex- 
pansion. Nevertheless, because the proposed potential is 
applicable for internuclear distances far beyond the equi- 



librium position (e.g., see Figs. 1 and 2), its asymptotic 
exponential behavior should not present an issue except 
for some extreme cases such as ultracold collisions. 

In conclusion, we have proposed an analytical three- 
parameter potential function for more than 200 weakly 
and strongly bound ground-state diatomics, including 
metastable molecular dications, with good accuracy over 
a significantly wide range of inter-nuclear distances. The 
determined potential energy curves and the associated 
vibrational levels agree well with literature data. When 
many-body effects are small, our simple pair potential 
function might also be useful in large-scale computer sim- 
ulations for complex systems. We anticipate that our 
model potential provides a useful guide towards supple- 
menting the potential curves obtained from the RKR and 
Dunham methods and a unified description of weakly and 
strongly bound diatomics. Extensions of our molecular 
orbital theory based approach to other types of diatomic 
systems are ongoing. 
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